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ABSTRACT 


H  prototype  coupled  Narine  Atmospheric  Boundary  Layer 
(NABL)  and  Oceanic  Boundary  Layer  (OBL)  model,  froa  physical 
nodels  which  have  been  separately  formulated,  is  developed. 
The  observational  efforts  will  be  directed  toward  coupling 
(local)  point  NABL  and  OBL  features.  Emphasis  in  the  char¬ 
acterization  of  local  features  will  be  on  the  evolutions  of 
the  adjacent  well-mixed  layers.  The  approach  will  be  to 
compare  observed  evolutions  in  the  oceanic  and  ataospheric 
boundary  layers  with  predictions  froa  bulk  models  wherein 
the  ocean  influences  the  atmosphere  through  the  surface  tem¬ 
perature  which,  in  turn,  is  influenced  by  atmospheric  forc¬ 
ing.  The  atmospheric  forcing  parameters  are  the  local 
surface  turbulent  fluxes  and  radiation.  Surface  turbulent 
fluxes  are  controlled  by  the  surface  temperature  and  atmos¬ 
pheric  mixed  layer  parameters.  Radiation  is  controlled  by 
predicted  cloud  formation  and  ii ssipatior..^ 
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I.  INTRODUCTION 


The  idea  that  both  atmospheric  and  subsurface  environ¬ 
mental  conditions  can  have  a  significant  effect  (enhancement 
or  degradation)  on  modern  weapon  and  sensor  systems  has  now 
become  a  reality  to  many  military  leaders.  This  has  led  to 
a  greater  demand  for  more  reliable  and  accurate  atmospheric 
and  oceanic  prediction  schemes.  Prediction  of  the  existence 
and  evolution  of  electromagnetic  ducts  are  important  in  pre¬ 
dicting  over  the  horizon  surface  to  surface  detection.  Pre¬ 
dictions  of  the  inversion  height  is  important  where  it  is 
known  that  optical  propagation  is  degraded  due  to  turbulence 
and  extinction.  Ocean  mixed  layer  predictions  are  important 
for  determining  convergence  zone  distances,  optimum  source 
level  depth,  and  cutoff  frequency  for  acoustic  surface  duct 
propagation. 

The  ability  to  provide  accurate  forecasts  requires  a 
complete  understanding  of  the  dynamical  processes  occurring 
in  each  of  the  atmospheric  and  oceanic  boundary  layers,  and 
as  important,  the  processes  associated  with  coupled  local 
changes  in  the  two  layers. 

Very  little  is  known  about  the  interactive  nature  of  the 
adjacent  turbulent  boundary  layers.  Th?  objective  of  this 
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thesis  is  to  formulate  a  model  which  couples  two  already 
existing  and  proven  models  for  the  respective  adjacent 
well-mixed  boundary  layers.  This  coupling  is  achieved  by 
matching  the  surface  boundary  conditions  for  stress,  mass 
flux  and  heat  flux  (Figure  1)  . 

This  coupled  model  may  then  be  used  to  examine  the  pro¬ 
cesses  of  coupled  local  oceanic  and  atmospheric  mixed  layers 
at  time  and  spatial  scales  of  12-18  hours  and  50-100  km, 
respectively.  Significant  changes  in  both  mixed  layers  can 
occur  at  these  scales.  Explaining  these  changes  requires 
coupled  evaluations  of  existing  atmospheric  and  oceanic 
well-mixed  layers  at  a  point. 

At  the  present  time,  we  only  have  qualitative  evidence 
*hat  changes  in  each  layer  are  correlated.  The  rigorous 
specifications  or  descriptions  of  the  coupling  will  require 
considerable  interpretive  efforts  and,  perhaps,  improved 


formulation  of  the  separate  models. 


II.  &&gtSfiOj£2 


A.  GENERAL  BOON  DART  LAYER  FEATORES 

The  Marine  Ataospheric  Boundary  Layer  (HABL)  extends 
froa  the  surface  through  the  capping  inversion.  The  Oceanic 
Boundary  Layer  (OBL)  extends  froa  the  surface  to  the  sea¬ 
sonal  theraocline.  The  air-sea  interface  is  bordered  by 
oceanic  and  ataospheric  turbulent  aixed  layers  which  effec¬ 
tively  insulate  the  quasi-geos  trophic  oceanic  and  atmos¬ 
pheric  flows  below  the  theraocline  and  above  the  inversion, 
respectively.  The  priaary  source  of  the  turbulence  is  the 
velocity  (current)  and  buoyancy  (density)  gradients  created 
by  interaction  a4-  the  interface.  The  large  vertical  mixing 
yields  nearly  homogeneous  (wsll-aixad)  wind,  temperature  and 
humidity  profiles  above  the  surface.  In  the  ocean,  the  term 
mixed  layer  refers  to  the  region  cf  vertical  uniformity  in 
mean  velocity  and  density.  Although  density  is  a  function 
cf  both  temperature  and  salinity,  it  is  usually  dependent 
upon  temperature  except  when  large  near-surface  salinity 
changes  occur  due  to  heavy  precipitation  or  ice  melt  in  high 
latitudes.  At  the  bottom  of  the  ocean  mixed  layer  and  at 


the  top  of  the  ataosphere  aired  layer  are  thin  transition 
regions  (thernocline  cr  inversion)  .  The  aired  layers  inter¬ 
act  with  the  tipper  and  lover  layers  at  these  regions  by 
means  of  turbulent  forced  entrainaent.  Typical  aodel  pro¬ 
files  depicting  the  boundary  layer  regions  for  both  the 
ataosphere  and  ocean  are  shown  in  Figure  2. 

The  well-mired  nature  of  both  the  turbulent  MABL  and  OBL 
implies  features  for  the  vertical  distribution  of  aean  val¬ 
ues  and  vertical  flures  of  wind,  temperature  and  humidity  in 
the  atmosphere  and  density  and  momentum  in  the  ocean.  One 
feature  is  that  the  properties  which  are  conserved  during 
mixing  can  be  treated  as  being  constant  with  height  (depth) 
in  the  boundary  layer.  A  second  feature  is  that  vertical 
flures  of  the  well-mixed  parameters  vary  linearly  with 
height.  These  features  enable  predictions  to  be  based 
solely  on  the  surface  and  inversion  (entrainment)  fluxes  if 
advection  is  negligible. 

Another  important  source  of  energy,  in  addition  to  tur¬ 
bulent  kinetic  energy,  in  both  mixed  layers  is  radiation, 
□alike  the  HABL,  most  of  the  solar  radiation  does  not  pene¬ 
trate  the  OBL.  Therefore,  downward  turbulent  heat  flux  in 
the  OBL  is  as  important  as  the  upward  flux.  Since  one  long 


15 


gure  2.  Simplified  atmospheric  aid  Oceanic  Boundary  Layer 
Tempera-cure  Profiles. 
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wave  (infrared)  radiation  flax  is  absorbed  and  emitted 
within  the  first  millimeter  of  the  sea  surface,  this  flux 
component  is  essentially  another  component  of  the  oceanic 
surface  heat  flux.  The  total  surface  heat  flux  is  comprised 
of  sensible  and  latent  heat  and  the  long  and  short  wave 
solar  radiation.  Beth  the  long  and  short  wave  (solar) 
fluxes  into  the  ocean  are  critically  dependent  upon  the 
existance  of  clouds  in  the  atmosphere. 

In  both  the  HABL  and  OBL,  buoyantly-created  velocity 
fluctuations  and  mixing  have  the  most  pronounced  effect  on 
mixed  layer  characteristics.  Garwood  (1977)  states  that 
buoyancy-forced  vertical  mixing  has  a  more  obvious  and 
direct  role  in  the  mechanical  energy  budget,  as  shown  in 
Figure  3.  In  the  atmosphere,  even  under  near  neutral  condi¬ 
tions,  buoyancy  forced  fluctuations  and  mixing  can  be  quite 
large.  These  buoyancy  driven  energetic  eddies  fill  the 
atmospheric  boundary  layer  from  the  surface  to  the  inversion 
where  they  entrain  warm,  dry  air  and  bring  more  momentum 
into  the  mixed  layer.  If  entrainment  causes  the  HABL  to 
extend  above  the  lifting  condensation  level,  clouds  or  fog 
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Figure  3.  Mechanical  Energy  Budget  for  the  Ocean  Mixed 
Layer. 

It  is  the  realization  that  buoyancy  driven  entrainment 
depends  on  the  near  surface  stability,  which  in  turn  is 
controlled  by  local  atmospheric  forcing  of  the  upper  ocean 
layer  which  makes  a  coupled  pra  diction  approach  necessary. 
A  circularity  in  cause  and  effect  arises  because  surface 
buoyancy  flux  induced  entrainment  nor  only  increases  the 
depth  of  the  MABL,  but  it  also  changes  its  effect  on  the 
ocean  mixed  layer.  The  role  of  the  surface  buoyancy  flux 
means  changes  ir.  the  oceanic  well-mixed  layer  could  be 


either  the  cause  or  the  result  of  entrainment.  Another 


cause  and  effect  relationship  is  the  formation  and  dissipa- 
tion  of  clouds  within  the  &BL.  Tlouds  can  be  caused  by 
changes  of  the  ocean  surface  tan perature.  In  turn,  clouds 
have  a  profound  effect  on  the  short  and  long  wave  radiation 
budget  for  the  OBL. 

These  cause  and  effect  relationships  shew  that  the 
assessment  of  relative  roles  of  dynamic  processes  in  both 
layers  requires  coupled  evaluations  using  relatively  coaplex 
physical  models.  1  data  set  collected  (Davidson  et  al  1983) 
over  a  48  hour  period  in  May  1978  demonstrates  the  niture  of 
coupled  local  changes  in  the  atmospheric  and  oceanic  bound¬ 
ary  layers.  For  a  period  from  1700  to  0500  local  time, 
clouds  existed  in  the  MABL  and  winds  increased  steadily  from 
6  to  9  as'1  during  this  period.  It  is  suggested  this 
increase  was  associated  with  deepening  of  the  MABL  (and 
therefore,  the  entrainment  of  additional  momentum  into  the 
mixed  layer)  which  is,  ir.  turn,  associated  with  the  surface 
temperature.  The  oceanic  boundary  layer  deepened  as  the 
wind  reached  a  maximum,  resulting  in  the  surface  temperature 
dropping  from  14.5°  to  13.0°C. 

There  were  other  atmospheric  response  and  forcing  fea¬ 
tures  observed  during  the  entire  48  hour  period  from 


5/19/1200  to  5/21/1200  PST  which  are  believed  to  have  been 
associated  with  coupling  between  the  two  layers.  Atmos¬ 
pheric  and  oceanic  sized  layer  observations  appear  in  Figure 
4. 

These  additional  features  are: 

a)  The  &BL  depth  increased  f ros  250  to  750  asters  over 
the  48  hour  period.  The  changes  occured  in  relatively 
short  intervals,  from  20/0000  to  20/0400,  and  froa 
21/0000  to  21/0500.  Tha  ABL  depth  reaained  at  near 
constant  depths  during  the  intervening  12  to  18  hour 
periods. 

b)  The  surface  layer  temperature  increase  from  20/0000 
to  20/1500  is  indicative  of  entrainment  of  overlying 
warm  air.  This  entrainment  presumably  was  a  factor 
in  the  wind  speed  increase  over  the  same  period.  The 
speed  increase,  in  turn,  presumably  caused  the  warm 
shallow  OBL  to  be  destroyed  which  dropped  the  surface 
temperature  after  20/1800  . 

c)  The  most  dramatic  change  ir.  the  3ABL  occured  during 
the  period  froa  21/0000  to  21/0500  which  was  immedi¬ 
ately  after  the  surface  temperature  had  dropped  due 


to  wind  mixing 


Questions  which  cculd  be  posed  on  the  basis  of  the  above 
data  are: 

a)  Did  the  drasatic  change  in  the  ABL  depth  at  21/0000 
occur  as  a  result  of  the  ocean  forcing  the  atmosphere? 

b)  Does  the  fact  that  the  aired  layer  evolved  with  rapid 
changes  followed  by  equilibrium  periods  indicate  feed* 
back  between  the  boundaries? 

c)  What  wculd  have  occurred  if  the  wind  had  reached  its 
aaxiaun  value  at  a  time  other  than  after  sunset? 

d)  Was  the  tine  of  occurrence  of  the  wind  maximum  and, 
hence,  the  ocean  aixing  controlled  by  the  coupled 
mixed  layer  system? 

This  thesis  will  attempt  to  answer  some  of  these  ques¬ 
tions  on  the  basis  of  the  formulated  coupled  model.  As 
stated  earlier,  all  answers  can  only  be  obtained  through 
considerable  interpretive  efforts. 

The  scope  of  the  overall  scientific  effort  is  illus¬ 
trated  in  Figure  5.  The  local  oceanic  descriptions  will  b» 
made  by  a  bulk  oceanic  mixed  layer  model  (the  Sarwood 
Model).  The  atmospheric  descriptions  will  be  mads  by  the 
Naval  Postgraduate  School  (NPS)  Marine  Atmospheric  Boundary 
Layer  model  (Davidson  et  al  198  3).  Tactically  significant 


operational  descriptions  (forecasts)  to  which  the  improved 
understandings  would  apply,  appear  on  the  extreme  right  hand 
side  of  Figure  5. 

B.  MARINE  ATMOSPHERIC  BOOMDARY  LAYER  (HABL)  MODEL 

The  UPS  boundary  layer  model  is  a  zero-order,  two  layer, 
integrated  aixed  layer  model.  The  two  layers  consist  of  the 
well-mixed,  turbulent  boundary  layer,  and  the  relatively 
non-turbulenc  free  atmosphere.  The  model  is  based  on 
entrainment  energetics  formulated  by  Stage  and  Businger 
(1981)  and  or.  radiative  transfers  described  by  Davidson  et 
al  (1983). 

The  inversion  (or  transition  zone)  separates  one  layer 
from  the  other.  In  a  zero-order  model,  the  inversion  is 
assumed  to  have  zero  thickness,  hence,  a  discontinuity,  or 
jump  occurs  ir.  the  profiles  of  the  conservative  parameters 
at  the  inversion. 

The  present  model  requires  as  input,  an  initial  atmos¬ 
pheric  sounding,  the  mean  winds  at  a  level  within  the  sur¬ 
face  layer  (10-30  meters)  and  the  surface  temperature. 
Boundary  layer  wind  values  for  the  duration  of  the  forecast 
period  can  be  obtained  from  larger  scale  models.  In  this 
study,  actual  hourly  wind  speed  measurements  were  input  as 
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part  of  the  model  initialization.  Sea-surface  temperature 
remains  unchanged  from  the  initial  input  value.  However,  in 
the  coupled  model,  the  sea-surface  temperature  will  change 
during  the  forecast  period  and  will  be  that  predicted  by  the 
OBL  model.  Hell-mixed  temperature  and  humidity  are  predic¬ 
ted  sc  the  surface  wird  and  wind  shear  at  the  inversion  are 
the  only  atmospheric  variables  which  have  to  be  prescribed. 
The  larger  scale  subsidence,  obtained  from  the  large  scale 
model  must  also  be  prescribed  for  the  forecast  period. 

Three  methods  car.  be  used  to  compute  the  large  scale 
vertical  velocity  (subsidence)  from  single  station  observa¬ 
tions.  These  three  methods  are  well  described  by  Gleason 
(1982)  and  include  the  kinematic  and  adiabatic  methods  and 
-he  integration  of  the  moisture  budget  equation  (Q-method)  . 
Gleason*  s  study  showed  the  Q-method  displayed  the  most  merit 
as  a  single-station  assessment  of  subsidence.  This  method 
was  used  to  compute  subsidence  for  this  study.  The  solar 
zenith  angle  used  in  the  computation  of  the  short  wave  radi¬ 
ative  flux  is  determined  from  the  latitude,  julian  day  and 
time  of  day  cf  the  initial  input. 

The  model  predicts  the  time  evolution  of  the  inversion 
height,  and  the  mixed  layer  values  of  3  (temperature)  and  g 
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(humidity),  as  well  as  their  jumps,  at  30-minute  intervals. 
The  prediction  model  is  run  out  to  24  hours.  The  steps  in 
the  prediction  computation  are  shown  in  Figure  6,  where  it 
is  ncted  that  procedures  are  the  same  for  clear  and  cloudy 
cases  except  for  entrainment  computation  and  estimation  of 
cloud  top  cooling. 

The  model  uses  the  standard  integrated  rate  equations 
(Tennekes  and  Dreidonks,  1981)  to  predict  the  time  evolution 
of  the  conservative  variables  and  their  jumps  at  the  inver¬ 
sion.  These  equations  are: 

h (Dx/Dt)  =  (w'xMq  -(w'x‘)h+  source  (1) 

h(Ddx/Dt)  =  hfx  (dh/di)  -  (w*  x*)0*  (w'xMh  •  source  (2) 

where  r  is  the  lapse  rate  above  the  inversion  and  the  source 
term  is  equal  to  -  (Fnh-  Fno)/pcp  for  x=9  and  equal  to  zero  for 
X=q.  is  the  net  radiative  flux.  The  subscripts  zero  and 

h  refer  to  surface  and  inversion  height  values,  respec¬ 
tively. 

Stage  and  Businger's  (1981)  entrainment  velocity  parame¬ 
terization  is  used  to  close  the  system  of  equations  and 
determine  the  time  evolution  of  the  inversion  height.  The 
closure  assumption  states  that  the  dissipation  rate  of  the 


igure  6.  Schematic  cf  Input,  Prescription  and  Computing 
Steps  in  H  ABL  Prediction. 
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turbulent  kinetic  energy  (TICE)  is  a  fraction  (1-A)  of  the 
preduction  rate.  '  A  1  is  the  entrainment  coefficient,  taken 
as  0.2. 

Mathematical  modeling  of  radiative  flux  transfer  is  an 
extremely  complicated  but  still  developed  subject.  Unfortu¬ 
nately,  even  the  simple  radiative  transfer  models  are  still 
extremely  complex  in  comparison  to  mast  other  physical 
parameterizations  used  in  the  mixed  layer  model.  Uncer¬ 
tainty  in  background  aerosols,  atmospheric  absorbing  gases 
(water  vapor,  carbon  dioxide,  ozone)  and  cloud  droplet  size 
spectra  are  sources  of  error  in  radiation  calculations. 
Long  and  short  wave  radiation  fluxes  are  computed  sepa¬ 
rately. 

The  long  wave  radiation  flux  calculation  was  modified  to 
permit  non-black  stratus  clouds  by  introducing  the  cloud 
emissivity,  £  ,  which  is  a  function  of  the  total  cloud  liq- 
uid  content,  W.  Since  the  cloud  liquid  water  profiles  are 
approximately  linear  with  height  (Davidson  et  al,  1982), 
cloud  liquid  water  content  and  cloud  emissivity  are  given 
by: 


0.5j?a(h  -  Zc)  q)h 

(3) 

£c  =  1  -  exp  (-aW) 

<“) 
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where  j?a  is  the  density  of  air  (1.25X  10'3  gm  cm*3  )  ,  h  is  the 
height  of  the  nixed  layer  (cloud  top),  Zc  is  the  lifting 
condensation  level  (cloud  botton) ,  a«0.158  (Slingo,  1982), 

and  g)h  is  the  liquid  water  content  at  the  cloud  top. 

The  long  wave  cloud  top  net  radiation  flux,  Lnh,  can  be 
calculated  from  the  cloud  top  temperature,  Th ,  and  the 
effective  radiative  sky  temperature,  Tsky  ,  using  the  Stefan- 
Boltzmann  Law.  At  the  cloud  bottom,  this  flux,  Lnc,  is  cal¬ 
culated  in  the  same  manner  using  the  sea  surface 
temperature,  Ts ,  and  cloud  bottom  temperature,  Tc .  These 
fluxes  are  given  by: 

Lnh  =ec<r(V-  Tsky4)  (5) 

Lnc*£c(T(Ts4-  Tc4)  (6) 

where  (T  is  Stefan’s  constant  (4.61  X  13"11  )  and  £c  ,  cloud 
emissivity,  is  obtained  from  equation  4.  The  net  long  wave 
flux  at  the  surface,  F|ong  ,  becomes: 

F.ong  =  Ws4-  <1-£c>Tsky4>  (7) 

where  T  is  the  average  temperature  of  the  cloud.  For  the 
cloud  free  case,  the  net  fluxes  are  calculated  at  Z  =  h  and 
Z=0  using  the  standard  method  (Fleagle  and  Businger,  1978) 
of  integrating  the  flux  emissivity  profile.  The  net  long 
wave  flux  at  the  surface  for  rhe  clear  sky  case  is: 
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P.ono  *  *u  *  ?d  <8> 

where  Fu  is  the  upward  radiative  flux  and  Fd  is  the  downward 
radiative  flux. 

Surface  flux  para  meterizatio  n  (bulk  aerodynamical)  for¬ 
mulae 


u*=  Cd  u10 
T*=  C e1/2(0o-  9) 
q*=  cj/2(q0-  q) 


are  used  to  determine  the  surface  fluxes  of  momentum,  sensi¬ 
ble  heat  and  latent  heat.  Th*sa  fluxes  are  given  by: 


2 

U  '  w*  a  u* 

(momen  turn) 

(10a) 

T '  w'  =  u*T* 

(sensible  ha  at) 

(10b) 

q»  w *  =  u *q* 

(latent  haat) 

(10c) 

where  Cd  and  C9  are  the  ten  mater  stability  dependant  drag 
coefficient,  9  is  the  potential  tamperatura,  and  q  is  the 
specific  humidity.  The  subscript  zero  denotes  the  surface 
value. 

The  shorn  wave  radiative  flux  is  calculated  using  the 
delta-Edingtcr  method  (Joseph,  1  976).  This  was  added  to  the 
model  to  account  for  heating  of  the  mixed  layer  by  solar 


33 


\*  -w 


1 '  l  u  .  1 

•  .  •  .  1  ■>- '  ' 


I 

radiation.  Short  nave  extinction  is  doeiQated  by  scattering 
(as  opposed  tc  absorption).  This  scattered  radiation  foras 
a  second  short  vave  radiative  component  usually  referred  to 
as  'diffuse'  solar  radiation.  Scattering  in  the  nixed  layer 
is  due  to  atmospheric  particles,  cloud  droplets  in  the 
cloudy  case  and  sea-salt  aerosols  in  the  cloud-free  case. 
This  MABL  model  computes  both  the  direct  and  diffuse  radia¬ 
tive  components  to  determine  the  total  short  vave  radiation 
flux  at  the  surface.  An  explanation  of  the  delta-Edington 
method  and  all  of  the  parameters,  atmospheric  factors,  and 
equations  involved  is  a  very  complex  subject  and  is  beyond 
the  scope  of  this  paper.  An  excellent  review  of  this  sub¬ 
ject  has  beer,  published  by  Pairall  (1981).  An  important 
short  wave  radiative  parameter  that  is  prescribed  in  the  ABL 
model  is  the  fraction  of  reflected  short  wave  radiation,  Ag, 
from  the  sea  surface.  It  is  taken  as  0.1. 

C.  OCEANIC  BOUNDARY  IAYEB  (OBL)  MODEL 

Garwood  (1977)  developed  an  ocean  mixed  layer  model 
using  the  Na  vier-Stokes  equation  of  motion  with  the  geostro- 
phic  component  eliminated,  the  continuity  equation  in  incom¬ 
pressible  water,  the  heat  equation  fron  the  first  law  of 
thermodynamics,  the  conservation  of  salt  equation,  and  a 


linearized  equation  of  state. 


Prediction  of  the  rate  of  3a opening  (or  retreat)  of  the 
mixed  layer  is  dependent  upon  an  understanding  of  the  dynam¬ 
ics  of  the  entrainment  process.  The  turbulence  of  the  over- 
lyinq  mixed  layer  provides  the  energy  needed  to  destabilize 
and  erode  the  underlying  stabLe  water  mass  (Garwood  1977). 
Therefore,  the  turbulent  kinetic  energy  budget  is  the  basis 
for  the  entrainment.  This  system  is  closed  using  a  mean- 
turbulent-field  modeling  of  the  vertically  integrated  equa¬ 
tions  for  the  individual  TKE  components,  plus  the  inclusion 
of  the  bulk  buoyancy  and  momentum  equations. 

Separate  vertical  and  horizontal  equations  for  TKE  are 
used  to  better  handle  the  mixing  process.  Buoyancy  produced 
energy  is  somewhat  more  efficient  than  shear  production  as  a 
source  of  energy  for  vertical  mixing  because  of  its  unique 
effect  on  the  vertical  component  of  rhe  turbulent  velocity. 
The  buoyancy  equation  is  generated  from  the  heat  and  salt 
equations  together  with  an  equation  of  state: 

}  =/o  -  90  )  +  $(S  -  S„  )]  (11) 
and  buoyancy  is  given  by: 

b=  q(  j>o  -/)//. 


(12) 


where  0  is  temperature,  S  is  salinity,  j>  density,  g  gravity, 
and  the  constants  «x  and  ^  are  the  expansion  coefficients  for 
heat  and  salt,  respectively.  The  *ilde  represents  the  total 
instantaneous  value  and  the  subscript  zero  denotes  an  arbi¬ 
trary,  but  representative,  constant  value.  Using  b  instead 
of  0  (tesperature)  for  the  definition  of  buoyancy  allows 
this  model  to  be  applicable  in  situations  where  evaporation 
and  precipitation  contribute  significantly  to  the  surface 
buoyancy  flux.  Temperature  is  in  most  cases  considered  to 
be  the  dominating  factor  in  the  density  profile.  The  rela¬ 
tive  importance  of  the  salinity  effect  on  the  short  term 
evolution  of  the  density  profile  is  generally  not  signifi¬ 
cant  (Miller,  1976). 

Temperature  and  salinity  profiles,  and  for  extended 
forecasts,  the  wind  driven  horizontal  current  profiles,  are 
required  as  input  for  model  initialization.  These  are  used 
to  compute  the  mixed  layer  depth.  The  Sarwocd  OBL  model 
defines  the  mixed  layer  depth,  h  ,  as  *:he  shallowest  depth 
at  which  the  observed  density  value,  a~[  ,  is  0.02  0^  units 
greater  than  the  observed  surface  density  value. 

Ocean  environmental  parameters  that  have  be  pre¬ 
scribed  include;  the  fraction  of  short  wave  radiation 


absorbed  in  the  top  one  seter  of  the  ocean,  the  radiation 
extinction  coefficient  for  determination  of  the  downward 
turbulent  heat  flux,  and  the  critical  Richardson  nuaber  for 
a  stability  adjustment  at  the  bottoa  of  the  aixed  layer. 

The  parameters  required  for  surface  boundary  condition 
computations  include  wind  speed  and  direction,  cloud  cover, 
sea  surface  temperature,  air  temperature  (dry  bulb),  dew¬ 
point  temperature,  incident  solar  radiation,  and  the  rate  of 
evaporation  (E)  and  precipitation  (P)  . 

The  turbulen^  fluxes  of  latent  heat,  Q0,  and  sensible 
heat,  Qh ,  are  estimated  using  the  bulk  aerodynamic  formu¬ 
lae  : 

Qe=  Cd  <.98ES  -  Ea)  u10  (13a) 

Qh=  Cd  (Ts  "  Ta)u10  (13b) 

and  the  net  back  radiation  is  estimated  from  the  empirical 
equation  (Husby  and  Seckel,  1  975  ): 

Qb=  1.14X  10*?(273.  16+TS)4  (.39  -  .5Ea/2)(1  -  .6C2)  (13c) 

where  Es  is  the  saturated  vapor  pressure  of  the  marine  air 
(0.98  corrects  for  salt  defects)  ,  Ea  is  the  vapor  pressure 
of  air  based  on  the  dew-point  temperature,  Ta  air  tempera¬ 
ture,  Ts  sea  surface  temperature,  and  C  is  the  fractional 
cloud  cover.  The  upward  heat  flux,  Qu,  is  then  given  as: 
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(14) 
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The  solar  radiation,  Q8,  is  given  by: 

Q9«  (1  -  aocb)  (1  -  .66C3)Q0  (15) 

where  Q0  is  the  clear  sky  radiation  given  by  (Seckel  and 
Beauday,  1973): 

Q0*  Aq*  A^cs  +  B1  sin  ^  ♦  i2cos  2$  *  B2sin  2 ^  (16) 

The  constants  a  and  b  are  adopted  from  Tabata  (1964)  and  the 
cubic  cloud  cover  correction  from  Laevastu  (i960).  The 

variable  is  the  aid-day  elevation  angle  of  the  sun.  The 
coefficients  (A0,At,  ect.)  were  calculated  by  harmonic  rep¬ 
resentation  of  the  values  predicted  in  the  Smithsonian 
Meteorological  Tables  (list  1958)  with 

*1=  (27r/365)  (t-  21)  (17) 

where  t  is  the  julian  day  of  the  year. 

Not  all  of  the  incoming  short  wave  radiation  penetrates 
the  ocean  mixed  layer.  Approximately  half  (for  open  ocean) 
is  absorbed  within  the  first  meter.  The  amount  absorbed 
varies  from  region  to  region  depending  on  the  amount  of 
absorbing  particles  such  as  phytoplankton  and  yellow  sub¬ 
stance.  More  radiation  will  be  absorbed  ir.  coastal  regions 
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than  in  the  open  ocean  regions.  This  portion  of  absorbed 
short  wave  radiation  is  therefore  considered  to  be  part  of 
the  upward  heat  flux.  The  reaainder  of  the  short  wave  radi¬ 
ation  does  penetrate  the  sixed  layer  and  is  attenuated  in  an 
exponential  fashion  depending  on  water  turbidity.  If  the 
fraction  of  solar  radiation  absorbed  in  the  first  meter,  Rf, 
is  given,  ther.  the  net  heat  at  the  surface  is  given  by: 

Qnet  *  Qu  ♦  (Hf)Q,  -  Qs  (18) 

From  these,  the  surface  fluxes  of  buoyancy  (heat  and 
salt)  and  momentum  can  then  be  computed.  The  turbulent  tem¬ 
perature,  salinity,  velocity  and  buoyancy  fluxes  are  given 
by : 


(T'  w*  ) 

=  "Qnat  ^fcp 

(19a) 

(S'  w' ) 

=  (P-E)  S0 

(19b) 

(u'  w*  ) 

=  U*2 

(19c) 

(b«  w* ) 

=  g[o<(T'w')  -  £(S'w')] 

(19d) 

respectively,  where  u*,  friction  velocity  in  air,  is  calcu¬ 
lated  by: 


"Cs  ~fa  Cd  U1  g 

(20) 

(21) 
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The  subscript  zero  refers  to  the  surface  value  and  XTS  is  the 

•  2 

surface  stress  (dynes  ca  ) .  In  this  lodel,  a  positive 
buoyancy  results  when  Qnet<0  and  E>P.  During  daytiae  peri¬ 
ods,  negative  buoyancy  flux  results  fro*  the  domination  of 
the  solar  heating  at  the  surface,  while  during  the  night 
periods,  a  positive  buoyancy  flux  is  produced  by  the  com¬ 
bined  effect  of  net  long  wavs  radiation  and  the  upward  tur¬ 
bulent  fluxes  of  heat  and  moisture. 

New  ocean  temperature,  salinity  and  wind  driven  current 
profiles  are  predicted  at  one  hour  intervals.  These  are 
then  used  to  predict  the  new  mixed  layer  depth.  The  steps 
in  the  prediction  computation  are  shown  in  Figure  7. 

D.  FORMULATED  COUPLED  BOUNDARY  LAYER  MODEL 

The  physical  processes  described  in  the  preceding  two 
sections  were  not  altered  during  the  coupling  steps.  Some 
modifications  were  involved  in  transferring  parameters 
between  the  models.  Also  necessary  was  the  uniformity  of 
units.  The  only  parameter  supplied  from  the  ocean  model  for 
use  in  atmospheric  computations  is  the  predicted  sea  surface 
temperature.  However,  this  prediction  is  influenced  by  many 


Initial  Profiles : 
Temperature ,  Sa 1 inity 


.gure  7 


Schematic  of  Inpur,  Prescription  and  Computing 
Steps  in  0 EL  Prediction. 
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Parameters  in  the  atmosphere  initialization  include; 
wind  speed  ar.d  direction,  dew-point  temperature,  latitude, 
julian  day,  and  time  of  day.  Short  wave  radiative  flux 
(direct  and  diffuse)  ,  long  wave  radiative  flux,  sensible  and 
latent  heat  fluxes,  and  the  surface  momentum  flux  were  com¬ 
puted  by  the  atmospheric  section  of  the  coupled  model  and 
input  to  the  ocean  model. 

The  first  problem  associated  with  the  coupling  involved 
‘he  different  prediction  time  steps  of  the  seperate  models. 
The  atmospheric  model  uses  a  30-minure  time  interval, 
whereas,  the  ocean  model  uses  a  ore-hour  time  step  interval. 
Assuming  no  significant  sea-surface  temperature  change  would 
take  place  ever  the  short  one-hour  interval,  the  ocean  mod¬ 
el's  needier- ion  scheme  was  inserted  into  the  atmospher® 
model  as  a  subroutine.  This  prediction  scheme  is  only 
called  every  c*her  30-minute  time  step.  The  ocean  model's 
initialization  routines  are  executed  immediately  upon  the 
completion  cf  the  atmosphere  model's  initialization. 

The  second  problem  was  with  wind  input.  The  A5L  model 
does  not  require  wind  direction  as  an  input  parameter,  only 
■‘•he  magnitude  of  the  wind  speed.  However,  the  ocean  model 
requires  wind  direction  tc  compute  the  horizontal  components 
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of  the  winds.  These  components  are  used  to  determine  the 

2 

ocean  turbulent  velocity  flux,  u*,  ,  components  used  in  the 
momentum  budget  computations,  k  new  subroutine  was  added  in 
the  initialization  scheme  to  determine  the  horizontal  wind 
components  from  wind  speed  and  direction  input.  These  com¬ 
ponents  are  determined  by: 

u10x  =  -  (sin  9)u10  (22a) 

uioy  =  “(cos  9)  u10  (22b) 

where  u10  is  the  specified  wind  speed  at  10  meters  and  the 
east-west  and  north- south  components  are  u10x  and  uioy  • 
respectively.  0  is  the  direction  from  which  the  winds  are 
blowing  relative  to  true  north.  Use  of  equations  20  and  21 
and  the  following  relationship: 


2 

allows  the  ocean  turbulent  velocity  flux  (u*w  )  tc  be  deter¬ 
mined  (assuming  the  density  of  seawater,  ,  to  be  approxi¬ 
mated  by  1)  by  first  computing  the  wind  stress  from 
atmospheric  parameters  and  then  equating  the  ocean  turbulent 

velocity  flux  to  the  computed  wind  stress  value.  The  hori- 

.  .  2 
zontai  components  of  the  turbulent  velocity  fluxes  u*wx  and 

2 

u*wy  ,  are  then  easily  computed  by: 


<2<») 

(25a) 

(25b) 


The  atmospheric  derived  values  are  computed  using  the  SKS 
unit  system  and  had  to  be  converted  to  the  CGS  system  prior 
no  any  oceanic  computations.  This  conversion  is  done  at  the 
start  of  each  ocean  prediction  run. 

As  discussed  in  section  B,  the  short  wave  radiative  flux 
is  calculated  using  the  delta-Ed ingt on  method.  The  long  wave 
radiative,  sensible  heat  and  latent  heat  fluxes  are  computed 
usinq  equations  10a,  10b  and  10c,  respectively.  A  unit  con¬ 
version  from  Watts  m  to  cal  cm  s  is  required  on  these 
values  prior  to  use  by  ocean  model.  This  conversion  is 
again  done  a*  the  start  of  the  ocean  prediction  run.  A  flow 
diagram  of  zhe  steps  in  the  coupled  model’s  prediction  com¬ 
putation  is  shown  in  Figure  8. 
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Figure  8.  Schematic  of  Atmospheric. and  Oceanic  Coupled 

Mode}..  The.  ocean  model  is  called  at  every  oth 
30-mmute  time  step  to  forecast  prescribed 
surface  temperature  for  the  MABL  prediction. 
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ill.  2A 21  iM2  E££2LIS 


The  data  sets  used  in  this  study  are  froa  the  Coopera¬ 
tive  Experiment  on  West  Coast  Oceanography  and  Meteorology 
from  21  September  to  12  October  1976  (CEBC3M-76)  and  froa  5 
to  23  May  1978  (CEWCOH-78)  .  Bath  were  conducted  in  the  Los 
Angeles-San  Diego  basin.  Oceanic  and  atmospheric  data  were 
collected  simultaneously  throughout  these  tiae  periods. 
XBT '  s  were  deployed  from  the  R/V  ACANIA  on  a  regular  basis 
(2-3  per  day)  in  conjunction  with  atmosphere  radiosondes. 
Sea  surface  temperature  and  atmospheric  data  such  as  winds, 
temperature,  and  relative  humidity  were  also  obtained  con¬ 
tinuously  through  these  periods  on  board  the  R/V  ACANIA. 

Three  24-hour  periods  were  selected  for  model  initiali¬ 
zation  and  verification.  The  first  two  cases  cover  a 
48-hour  period,  5/19/0500  to  5/21/0500  PST,  during 
CEWCOM-78.  This  pericd  coincides  with  the  observations  dis¬ 
cussed  in  chapter  2,  in  which  coupled  air-sea  interactions 
were  cited  as  possible  explanations  for  changes  in  both  the 
atmospheric  and  oceanic  mixed  layer.  Satellite  imagery 
showed  uniformly  increasing  stratus  coverage  (thin  to  heavy) 
during  rhis  48-hour  period.  A  clear  sky  situation  was 
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chosen  for  the  third  case.  This  24-hour  period,  10/04/1900 
to  10/05/1900  PST,  was  during  CBWCOH-76. 

The  general  location  of  the  R/V  ACAHIA  during  these 
periods  is  shown  in  Figures  9,  10  and  11.  During  case  1  the 
ship  was  at  anchor  for  the  majority  of  the  time  period.  The 
anchorage  location  is  shown  in  Pigure  9.  For  cases  2  and  3 
the  R/V  ACAHIA  was  cruising  slowly  (2-3  knots)  returning  to 
an  initial  point  approximately  every  12  hours.  winds  were 
primarily  from  the  west  during  the  case  1  period.  The 
winds  shifted  slightly  to  a  west-northwest  direction  midway 
through  the  case  2  period.  Winds  for  case  3  were  primarily 
from  the  northwest  and  were  much  slower  than  in  either  case 
1  or  2. 

Two  predictions  were  analyzed  for  each  case,  one  predic¬ 
tion  was  based  or.  the  coupled  air-sea  boundary  layer  model, 
and  the  second  prediction  was  with  the  atmospheric  model 
with  a  fixed  surface  temperature.  A  separate  oceanic  pre¬ 
diction  using  fixed  atmospheric  parameters  (persistence)  was 
not  made.  This  was  the  result  of  an  initial  modification 
that  removed  the  wind  stress  and  heat  flux  computation  pack¬ 
age  in  the  ocean  model  in  preparation  for  coupling  with  the 
atmosphere  model. 
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Figure  9.  Anchorage  Position  of  P/7  ACA5frA,ir.  the  7icini*y 
of  San  Hicolas  Island  from  5/19/7S  to  5/2o/^8.  1 


Presentation  of  results  will  be  in  two  parts.  The  first 
ro  be  presented  is  the  coupled  model's  performance  with 
respect  tc  observed  oceanic  and  atmospheric  variations.  The 
second  will  deal  wi*h  a  comparison  of  predicted  values 
between  the  coupled  model  and  the  atmospheric  model.  All 
-abies  and  figures  are  grouped  together  by  case  number  a*, 
the  end  of  the  chapter. 


« 


5 


Pt. 

Conception 


CEWCOM-78 


Figure  10.  Positions  of  H/V  ACANIA  from  5/18/78  to  5/25/^8 
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A.  RESULTS  FROM  COUPIED  MODEL 

1.  C£5E  J  (0550  15  Maz  to  2522  2D 

A  listing  of  predicted  atmospheric  and  oceanic  aired 
layer  quantities  for  this  case  appear  in  Table  I.  Values 
are  for  each  30-minute  time  step  and  since  oceanic  values 
are  only  computed  at  one-hour  intervals,  the  same  value 
appears  at  two  time  steps.  The  parameters  listed  include: 
height  of  inversion  (Zi)  ,  height  of  lifting  condensation 
level  (Zlcl)  ,  wind  speed  (Hind),  atmospheric  mixed  layer 
potential  temperature  <Th) ,  evaporation  duct  height  (Ze)  , 
sea  surface  temperature  (SST)  ,  ocean  mixed  layer  depth 
(MLD) ,  and  the  net  ocean  surface  heat  flux  (TW)  . 

Predicted  values  are  also  shown  in  Figures  12 
through  17.  Figure  12  shows  the  time  evolution  of  predic¬ 
ted  sea-surface  temperatures  (solid  line)  and  the  values  of 
the  measured  surface  temperatures  (dashed  line).  Figure  13 
shews  the  time  evolution  of  the  ocean  mixed  layer  depth. 
Fiqures  14-16  show  ocean  vertical  profiles  for  temperature, 
u  and  v  wind  driven  current  velocities,  and  the  salinity 
change  from  an  initial  value  of  33.5  ppt,  respectively. 
Figure  17  shows  predictions  of  atmospheric  parameters. 
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High  solar  insulation  luring  the  daytiae  results  in 
a  net  downward  heat  flux  (ocean  warning)  .  This  can  be  read¬ 
ily  observed  in  the  net  surface  heat  flux  values  (TW)  ,  in 
Table  I,  with  the  negative  sign  appearing  during  the  daytiae 
hours.  Maximum  heating  of  the  ocean  occurred  between  1300 
and  1400  local  tiae.  This  saae  pattern,  warning  in  the  day¬ 
tiae  hours,  is  also  observed  in  the  sea  surface  teaperature 
values.  Table  I  and  Figure  12. 

The  winds  appeared  to  be  “airly  constant,  8ns'1, 
between  0500  and  1100  on  the  nineteenth,  find  nixing  coupled 
with  convective  cooling  at  the  surface  caused  the  mixed 
layer  depth  to  deepen  slightly  froa  19.5  to  20.9  m  (Figure 
13).  This  deepening  and  net  heat  loss  at  the  surface  caused 
the  surface  teaperature  to  drop.  As  heating  of  the  sea  sur¬ 
face  began  to  increase  due  to  the  increase  of  solar  radia¬ 
tion,  a  shallow  mixed  layer  of  14  meters  formed. 

The  winds  rapidly  decreased  from  9  to  5  ms’1  at  1200. 
This  coincided  with  a  near  maximum  value  of  radiative  heat¬ 
ing  causing  a  shallow  7  meter  mixed  layer  to  form  in  less 
than  1  hour.  As  daytime  heating  continued,  a  shallower 
mixed  layer  cf  3.5  meters  formed.  This  process  is  referred 
to  as  the  'afternoon  effect'  in  acoustic  literature  and  is  a 
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quite  common  phenomenon  in  many  areas  of  the  ocean.  In  Fig¬ 
ure  14,  the  7  meter  mixed  layer  is  barely  discernible  in  the 
1300  vertical  temperature  profiLe.  The  weak  gradient  at  the 
base  of  the  mixed  layer  only  required  weak  mixed  layer  ero¬ 
sion  (increase  in  winds  or  convective  cooling)  to  break  down 
the  shallow  mixed  layer.  During  the  late  afternoon,  the 
winds  began  *0  increase,  and  with  the  decreasing  of  solar 
radiation,  the  shallow  thermocline  was  eroded  away.  The  MLD 
continued  deepening  for  the  remainder  of  the  time  period  due 
to  wind- driven  and  convective  mixing. 

Since  sea-surface  temperature  is  the  sole  input  from 
the  ocean  model  to  the  atmosphere  model,  a  comparison  is 
made  of  predicted  and  actual  sea-surface  temperature  values 
(Fiqure  12).  Since  the  R/V  ACANTA  was  at  anchor  during  the 
majority  of  this  time  period,  this  is  considered  the  best 
case  to  determine  ho*  well  sea  surface  temperature  changes 
due  to  atmospheric  forcing  are  being  predicted.  The  coupled 
model  predic*ion  yielded  a  3.25 °C  temperature  increase  dur¬ 
ing  the  daytime  heating  period  with  a  maximum  temperature  of 
15 . 4°C  at  1600  local  time.  A  maximum  surface  temperature  of 
15.6°C  was  observed  at  1600  giving  a  0.2°C  error. 


After  the  daytime  heating  period,  the  surface  tem- 
peraturas  began  to  drop.  The  ainiaua  predicted  surface  tem¬ 
perature  at  the  end  cf  the  24-hour  period  was  15.2°C  which 
yields  a  0. 2°C  teaperature  drop.  The  0.06°C  temperature 
drop  prior  to  ocean  warming  lakes  the  magnitude  of  the  pre¬ 
dicted  diurnal  warming  and  cooling  to  be  approximately 
equal.  The  actual  minimum  temperature  at  the  end  of  the 
period  was  14.96°C.  Although  there  are  small  errors  in  the 
magnitude  of  the  diurnal  heating  and  cooling  affects,  the 
time  evolution  of  this  phenomenon  is  predicted  quite  well. 

Precipitation  was  not  observed  during  this  period, 
therefore  the  surface  moisture  flux  and  surface  salinity 
flux  were  due  entirely  to  evaporation.  The  slight  predicted 
increase  in  salinity  from  33.5  to  33.5038  ppt,  shewn  in  Fig¬ 
ure  16,  is  due  to  evaporation. 

An  unstable  condition  existed  in  the  atmosphere  dur¬ 
ing  this  prediction  period.  A  strong  buoyancy  flux  at  the 
surface  increases  entrainment  and  should  cause  the  inversion 
height  to  increase  unless  subsidence  (associated  with  shal¬ 
lowing  of  the  inversion  height)  is  strong  enough  to  prevent 
this  increase.  In  this  case,  subsidence  was  relatively  weak 
(-.0028  ms"1)  and  entrainment  was  able  to  increase  the  inver¬ 
sion  height. 
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Cloud  occurance/dissipat  ion  was  of  interest  because 
cf  its  affect  on  radiation  at  the  surface.  With  the  lifting 
condensation  level,  LCL,  lower  than  the  inversion  throughout 
this  time  period,  stratus  clouds  should  be  present  during 
the  entire  period.  This  feature  was  observed.  During  this 
period  the  LCL  decreased  from  2  43  to  125  aeters  while  the 
inversion  height  increased  from  256  to  348  meters.  This 
would  indicate  a  thickening  of  the  clouds  over  this  24-hour 
period  which  was  also  observed  to  be  true.  The  air-sea 
coupling  effect  on  predicting  these  cloud  features  cannot 
be  identified  until  the  results  with  a  constant  sea  surface 
temperature  forcing  are  examined.  However,  on  the  basis  of 
the  predicted  temperature  change,  a  significant  influence  by 
the  coupling  is  not  expected. 

2.  CASE  2  (0  500  20  May  to  0  500  21  May) 

A  listing  of  predicted  atmospheric  and  oceanic  mixed 
layer  quantities  for  this  period  are  given  in  Table  III  and 
shown  graphically  in  Figures  19  through  24.  As  in  case  1, 
daytime  heating  at  the  surface  is  apparent  by  the  negative 
siqns  in  front  of  the  net  surface  heat  flux  values.  The 
maximum  downward  heat  flux  value  was  not  as  high  as  in  case 
1  and  was  probably  due  to  thicker  cloud  conditions  for  this 
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case.  However,  increasing  surface  temperatures  were 
observed  during  this  heating  period. 

The  winds  increased  steadily  from  6  to  11  ms"1 
throughout  the  first  half  of  this  psriod  and  slowly 
decreased  to  about  9.4  ms'1  at  the  end  of  the  period.  The 
higher  winds  should  be  associated  with  increased  wind  mixing 
cf  the  ocean  mixed  layer.  The  increased  winds  caused  the 
mixed  layer  depth  to  deepen  at  the  end  of  this  peiod.  This 
also  reduced  the  shallowing  of  the  mixed  layer  during  the 
afternoon  radiative  heating  period.  The  shallowest  mixed 
layer  formed  at  16  meters  versus  3.5  maters  in  the  previous 
case.  In  Figure  21,  the  20  mater  thermorline  is  seen  tc  be 
very  weak,  and  with  the  higher  winds,  tha  shallower  after¬ 
noon  mixed  layer  lasted  only  a  couple  of  hours. 

The  effects  of  the  diurnal  heating  and  cooling  peri¬ 
ods  caused  a  0.1°C  temperature  increase  and  a  0.15°C  temper¬ 
ature  decrease,  respectively.  Both  valuas  ara  also  smaller 
than  those  predicted  in  case  1 .  Reasons  for  the  smaller 
predicted  increase  are:  1)  thicker  clouds  reduced  tha  amount 
cf  solar  heading  during  the  day,  tharby  reducing  the  diurnal 
temperature  increase;  2)  the  higher  winds  increased  the  ver¬ 
tical  mixing,  causing  cooler  water  below  to  be  mixed  into 


the  upper  layer  and  also  increased  the  latent  heat  loss  due 
to  evaporation,  thus  reducing  the  total  net  heat  flux  at  the 
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surface;  and,  3)  the  larger  air-sea  temperature  difference 
increased  the  sensible  heat  loss,  which  also  reduced  the  net 
heat  flux.  The  decrease  in  the  downward  solar  radiative 
heating  and  increase  in  sensible  and  latent  heat  losses  were 
somewhat  balanced  by  the  decrease  in  the  long  wave  heat  loss 
due  to  the  presence  of  the  cLouds.  This  decrease  in  the 
long  wave  heat  loss  was  also  the  primary  reason  for  the 
slightly  smaller  value  in  the  nighttime  cooling  period. 

Comparisons  cf  the  predicted  surface  temperature 
values  with  the  actual  values  (Figure  19),  show  a  noticeable 
difference.  This  is  probably  due  to  the  fact  that  at  0700 
local  time  or.  20  May,  the  R/7  ACANIA  began  cruising  slowly 
to  the  northwest.  The  regional  horizontal  temperature 
structure  must  be  taken  into  account.  The  southern  Califor¬ 
nia  coastal  bathymetry  and  associated  currents  greatly 
influence  the  sea  surface  temperature  structure.  The  Cali¬ 
fornia  Current  brings  in  colder  water  from  the  north  as  it 
heads  south  along  the  coast.  Near  Point  Conception,  the 
current  begins  to  turn  to  the  southwest  a*ay  from  the  coast. 
Due  to  coas+al  bathymetry.  Point  Conception  also  provides 


sheltering  for  the  intrusion  of  warm  water  along  the  coast 
from  th9  south  by  the  Davidson  Current.  The  eddies  and 
meanders  of  cold  and  warm  water  make  the  sea  surface  temper¬ 
ature  structure  complex  in  the  CEWCOM  experiment  areas.  The 
surface  temperature  variation  observed  by  the  B/V  ACANIA 
appears  to  be  due  to  oceanographic  features. 

The  1°C  temperature  increase  in  the  first  two  hours, 
when  the  ship  began  moving  to  the  northwest,  indicates  that 
the  ship's  ■‘■ransit  was  through  one  of  the  small  warm  eddy 
features.  Two  reasons  for  this  are  1)  Radiative  heating  of 
xhe  ocean  did  not  even  begin  until  after  the  observed  sur¬ 
face  temperature  had  already  begun  to  rapidly  decrease,  and 
2)  the  higher  winds  would  tend  to  mix  the  surface  waters 
with  the  cooler  water  from  below,  thus  causing  a  decrease  in 
temperature  not  an  increase.  The  short  time  duration  of 
this  warm  feature  also  suggests  a  transit  through  an  ano¬ 
maly.  With  the  ship  only  cruising  at  2-3  knots,  the  warm 
anomaly  still  only  lasted  less  than  4  hours. 

The  ship  turned  on  a  northerly  course  at  1300  and 
-.he  rapid  temperature  decrease  stopped.  It  is  reasonably 
safe  to  infer  from  the  sudden  change  in  the  surface  tempera¬ 
ture  trend  (rapid  decrease  to  slight  increase)  that  the  ship 
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is  initially  crossing  the  temperature  gradient  contours  and 
then  followed  a  parallel  course  with  these  contours.  This 
would  make  oceanic  conditions  more  compatible  with  those  in 
case  1.  Since  the  0.4°c  temperature  increase  daring  the 
daytime  heating  period  compares  favorably  with  the  value 
observed  while  the  ship  was  at  anchor,  the  increase  is 
therefore  presumably  due  to  the  diurnal  heating. 

The  ship  turned  back  to  a  west-northwest  course 
shortly  after  1800  and  the  surface  temperature  immediately 
began  to  decrease  1.5°C  in  4  hours.  Two  hours  after  the  end 
of  this  period,  0700,  the  ship  came  about  and  paralleled  its 
previous  course  but  in  the  opposite  direction.  An  examina¬ 
tion  of  the  observed  sea  surface  temperatures  shows  a  rapid 
increase  comparable  to  the  decrease  while  heading  in  the 
opposite  direction.  This  further  substantiates  the  observed 
large  temperature  changes  are  not  due  to  atmospheric  forc¬ 
ing,  but  rather  to  the  ship's  movement  through  an  oceano¬ 
graphic  eddy  or  front. 

Unstable  atmospheric  conditions  existed  during  this 
period  as  well.  The  1.5°C  drop  in  the  air  temperature  and  a 
much  smaller  drop  in  the  ocean  temperature  (0.2°C)  coupled 
with  higher  wind  speeds  resulted  in  a  large  surface  buoyancy 


flux.  This  resulted  in  a  much  higher  predicted  inversion 
height  change  due  to  entrainment.  The  slightly  smaller  sub¬ 
sidence  rate  was  also  a  factor  in  the  greater  predicted 
increase  of  the  inversion  height. 

The  lifting  condensation  level  fluctuates  up  and 
down  unlike  the  steady  decrease  in  height  in  case  1.  How¬ 
ever,  the  height  of  the  LCL  never  increased  above  the  ini¬ 
tial  height  and  so  with  the  large  increase  in  height  of  the 
inversion,  thickening  of  the  stratus  cloud  deck  is  predicted 
to  continue  during  this  period  as  wall. 

3.  CASE  3  (JSOO  oa  OCT  to  1900  05  OCT) 

A  listing  of  predicted  atmospheric  and  oceanic  mixed 
lay^r  quantities  are  given  in  Table  V  and  shown  graphically 
in  Figures  26  through  31.  Although  the  time  period  for  this 
case  is  in  October,  there  are  still  many  features  occuring 
that  are  similar  to  the  first  two  cases  in  May.  The  primary 
differences  in  this  case,  other  than  being  a  cloud-free 
case,  are  1)  the  air  and  surface  temperature  are  much  warmer 
and  2)  a  much  higher  value  of  subsidence  was  observed. 

Daytime  heating  is  apparent  in  predictions  for  the 
clear  sky  case,  as  expected.  The  magnitude  of  •‘■he  net  ocean 
surface  heat  flux  is,  however,  smaller  than  the  value  for 
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the  thin  stratus  deck  condition  in  case  1.  This  is  presun* 
ably  due  to  the  difference  in  seasonal  heating  rate  and  will 
not  be  considered  part  of  the  local  coupling  effects. 

The  winds  in  case  3  ace  such  weaker  than  in  the 
first  two  cases.  with  the  light  winds  and  high  downward 
surface  heat  flux  during  the  day,  a  nuch  shallower  nixed 
layer  might  be  expected.  The  ocean  mixed  layer  depth  (HID) 
values  in  Table  V  (also  shown  graphically  in  Figure  25)  show 
that  an  ocean  nixed  layer  of  1  meter  suddenly  forms  in  less 
*han  an  hour  after  heating  began  at  0930.  This  shallow 
thermoclire  is  much  stronger  than  in  either  case  1  or  2,  and 
persists  throughout  most  of  the  day  (Figure  28).  As  the 
winds  slowly  increased  late  in  the  afternoon,  the  shallow 
mixed  layer  deepens.  Without  winds  to  mix  this  shallow 
region,  surface  heating  is  much  greater  in  this  case.  A 
0.6°C  surface  temperature  increase  was  predicted  to  occur. 

As  in  case  2,  the  ship’s  movement  played  a  major 
role  in  the  large  variations  in  observed  surface  temperature 
values.  This  prevented  the  verification  of  the  predicted 
surface  temperature  increase  during  the  daytime  heating 
period.  Figure  11  shows  the  complex  horizontal  temperature 
structure  in  the  CEWCOM-76  experiment  area  for  the  first 


week  of  October.  Fro*  this  figure,  it  is  hard  to  explain 
the  large  2.8°C  teaperature  drop  during  the  first  half  of 
th9  period  as  the  ship  appeared  to  be  heading  into  warmer 
water.  The  rapid  temperature  increase,  starting  at  0900,  is 
obviously  due  to  the  ship  transiting  into  the  war*  eddy  fea¬ 
ture  north  of  Santa  Catalina  Island  (Figure  11). 

Unstable  atmospheric  conditions  are  also  present 
during  this  time  period.  The  high  subsidence  caused  the 
inversion  to  be  forced  downward  to  within  72  meters  of  the 
surface.  The  lifting  condensation  level  decreases  even  more 
rapidly  and  finally  becomes  lower  than  the  inversion  at  the 
last  hour  of  the  forecast  period.  This  implies  the  forma¬ 
tion  of  a  low  stratus  cloud  deck.  From  Figure  33,  fog  and 
stratus  clouds  were  observed  to  form  in  the  CEWCOM  experi¬ 
ment  area. 

B.  COMPARISON  OF  COUPLED  AND  ATMOSPHERIC  MODEL 

The  second  and  final  phase  of  this  study  is  to  compare 
the  previous  coupled  air-sea  boundary  layer  model  results 
with  the  atmospheric  model  results  using  a  constant  sea  sur¬ 
face  teaperature.  This  comparison  is  to  determine  if  any 
differences  existed  between  the  atmospheric  predictions,  and 
if  so,  determine  if  these  differences  are  due  to  coupled 
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interactions.  A  listing  of  the  atmospheric  model  predic¬ 
tions  of  atmosphere  mixed  layer  parameters  for  cases  1,  2 
and  3  are  found  in  Tables  II,  17  and  71,  respectively. 
Graphical  presentations  of  these  results  are  shown  in  Fig¬ 
ures  18,  25  and  32.  All  tables  and  figures  are  grouped 
together  by  case  number  at  the  end  of  this  chapter. 

The  initial  comparison  reveals  no  significant  difference 
in  the  two  predicted  sets  of  mixed  layer  parameters,  indi¬ 
cating  a  negligible  effect  of  coupled  interactions  on  atmos¬ 
pheric  predictions  for  this  short  time  period.  The 
predicted  inversion  heights  are  essentially  identical  in 
both  sets  of  predictions  for  all  three  cases  (maximum  dif¬ 
ference  of  2  meters)  .  This  is  also  found  to  be  true  for 
humidity  and  even  the  evaporation  duct  height  predictions. 
This  observation  is  somewhat  surprising  since  there  are  sig¬ 
nificant  interactions  between  the  surface  fluxes  (evapora¬ 
tion)  and  these  parameters.  The  close  proximity  of  these 
parameters  to  the  ocean  surface  would  also  imply  a  more  sig¬ 
nificant  impact  even  with  small  changes  in  sea-surface  cem- 
pe  ratures. 

The  only  parameter  that  seems  to  be  clearly  affected  by 
the  coupling  is  the  lifting  condensation  level.  However, 


s 

I 

I 

1 

I 


even  for  the  cloudy  conditions  of  cases  1  ar.d  2,  this  shows 
only  a  slight  deviation.  A  maximum  difference  of  only  5 
meters  is  observed  between  the  predictions  for  these  cases. 
In  the  clear  sky  case  under  warmer  and  tore  huaii  condi¬ 
tions,  there  is  a  30  meter  difference  observed  in  the  two 
LCL  predictions.  However  this  occurred  as  the  atmosphere 
approach  near  neutral  conditions.  This  may  indicate  the 
model's  inability  to  make  reliable  forecasts  under  stable 
atmospheric  conditions. 

Since  a  comparison  of  predicted  ocean  parameters  between 
a  coupled  ar.d  independent  ocean  model  was  net  made,  no  con¬ 
clusions  on  this  topic  can  be  made  at  this  time.  Although 
tentativa  results  indicate  no  definite  advantages  using  a 
coupled  model  for  the  atmospheric  predictions,  the  reverse 
may  very  well  be  true  for  oceanic  predictions. 
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IV.  COWCmSIONS/R  ECOWMSNDMIOHS 


On  the  basis  of  a  limited  study  of  several  24-hour  pre¬ 
dictions,  surface  temperature  variations  are  not  large 
enough  to  require  ocean  mixed  layer  processes  in  short  term 
atmospheric  forecasts.  Carefully  selected  data  sets  and 
considerable  interpretive  efforts  are  essential  if  general 
conclusions  are  to  be  drawn  on  the  relative  roles  of  dynamic 
processes  in  bo*h  boundary  layers.  Further  studies  are 
required  to  determine  the  model's  ability  to  better  handle 
forecasts  at  longer  time  scales.  Further  studies  should 
also  include  sensitivity  analyses  to  define  situations  for 
which  air-sea  interaction  effects  are  greatest.  Studies  of 
oceanic  sensitivity  to  atmospheric  forcing  are  required  to 
determine  the  impact  it  will  have  on  acoustical  forecasts. 

The  CEWCOM-76  and  78  data,  although  able  to  supply 
simultaneous  and  oceanic  data  ,  suffered  in  ‘wo  major 
respects  First,  the  complexity  of  the  horizontal  temperature 
structure  with  meanders  of  warm  and  cold  eddy  features  had 
definite  impacts  on  observed  data.  Second,  data  was  col¬ 
lected  while  the  R/Y  ACANIA  was  transiting  through  a  complex 
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ocean  temperature  structure.  Changas  in  observed  conditions 
under  these  circumstances  can  and  did  mask  the  interaction 
processes  of  interest. 

Although  the  primary  results  of  this  study  indicate 
coupled  air-sea  models  are  probably  not  necessary  for  short 
term  atmospheric  forecasts,  a  coupled  model  still  has  advan¬ 
tages,  especially  if  long  ten  forecasts  are  to  become  pos¬ 
sible.  Such  a  model  would  not  be  the  single  one-dimensional 
one  considered  here.  However,  this  model  offers  many  poten¬ 
tial  benefits  as  a  multi-purpose  model  with  single  station 
assessment  and  prediction  capabilities  in  the  absence  of 
shore-based  support.  A  sound  velocity  profile  (SVP)  predic¬ 
tion  capability  could  be  easily  incorporated  into  this 
model.  In  conjunction  with  present  prediction  capabilities, 
a  multi-purpose  model  would  be  available  for  a  wide  range  of 
needs.  Ocean  mixed  layer  depth  and  SVP  predictions  can  be 
used  in  predicting  ASM  sensor  capabilities.  Environmental 
affects  cn  the  performance  of  weapons  and  sensors  that 
depend  or.  electromagnetic  and  electro-optical  wave  propaga¬ 
tion  can  be  evaluated  using  the  same  model. 
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TABLE  I 


Coupled  Model  Predicted  A*aospheric  and  Oceanic  Mixed  Lave: 

Quantities-  Case  1. 


HOUR 

Zi 

Zlcl 

WIND 

Th 

SST 

MLD 

TU 

Ze 

0430 

256 

243 

0.0 

13.3 

15.2 

0.0 

0 . 0000 

0.0 

0500 

252 

216 

8.0 

13.4 

15.2 

19.5 

0.0000 

7.2 

0530 

250 

174 

8.0 

13.3 

15.2 

20.1 

.  0067 

6.9 

0600 

250 

130 

8.0 

13.0 

15.2 

20.1 

.0067 

6.4 

0630 

252 

116 

8.0 

12.9 

15.1 

20.4 

.0  028 

6.0 

0700 

256 

121 

8.1 

12.8 

15.1 

20.4 

.0028 

6.0 

0730 

259 

127 

8.1 

12.9 

15.1 

20.7 

.0011 

6.1 

0800 

262 

132 

8.1 

12.9 

15.1 

20.7 

.0011 

6.1 

0830 

265 

136 

8.1 

12.9 

15.1 

20.3 

-.0007 

6.2 

0900 

267 

139 

8.1 

12.9 

15.1 

20.8 

-.0007 

6.2 

0930 

270 

142 

8.1 

12.9 

15.2 

20.9 

-.0028 

6.3 

iOOO 

272 

145 

8.2 

13.0 

15.2 

20.9 

-.0028 

6.3 

1030 

274 

148 

8.2 

13.0 

15.2 

16.4 

-.0045 

6.3 

1100 

276 

151 

8.2 

13.0 

15.2 

16.4 

-.0045 

6.3 

1130 

278 

153 

6.6 

13.1 

15.2 

14.  1 

-.0057 

6.3 

1200 

279 

153 

5.0 

13.0 

15.2 

14.1 

-.0057 

5.7 

1230 

280 

153 

5.1 

13.0 

15.3 

7.2 

-.0066 

5.2 

1300 

281 

153 

5.3 

12.9 

15.3 

7.2 

-.0066 

5.2 

1330 

282 

153 

5.4 

12.8 

15.3 

3.6 

-.0066 

5.2 

1400 

283 

153 

5.6 

12.8 

15.3 

3.6 

-.0066 

5.3 

1430 

284 

152 

5.8 

12.8 

15.4 

5.  1 

-.0052 

5.4 

1500 

285 

150 

5.9 

12.7 

15.4 

5.1 

-.0052 

5.4 

1530 

286 

143 

6.1 

12.7 

15.4 

6.5 

-.0032 

5.5 

1600 

288 

145 

6.3 

12.6 

15.4 

6.5 

-.0032 

5.6 

1630 

289 

142 

7.0 

12.6 

15.4 

3.3 

-.0008 

5.6 

170  0 

291 

138 

7.8 

12.5 

15.4 

8.3 

-.0008 

5.9 

1730 

293 

134 

7.7 

12.5 

15.3 

12. 0 

.0014 

6.1 

1800 

296 

130 

7.7 

12.5 

15.3 

12.0 

.0014 

6.1 

1830 

299 

127 

7.6 

12.4 

15.3 

14.8 

.  0029 

6.1 

1900 

301 

125 

7.6 

12.4 

15.3 

14.8 

.0029 

6.1 

1930 

304 

125 

7.5 

12.4 

15.3 

17. 0 

.0034 

6. 1 

2000 

307 

125 

7.9 

12.4 

15.3 

17.0 

.  0034 

6.1 

2030 

310 

127 

8.3 

12.4 

15.3 

18.7 

.  0034 

6.2 

2100 

313 

128 

8.1 

12.4 

15.3 

18.7 

.0034 

6.4 

2130 

316 

i29 

7.9 

12.3 

15.3 

20 . 6 

.  0033 

6.3 

2200 

319 

129 

7.7 

12 . 3 

15.3 

20.6 

.0038 

6.3 

2230 

322 

129 

7.5 

12.3 

15.3 

21 . 1 

.  0036 

6.2 

2300 

324 

i29 

7.3 

12.3 

15.3 

21 . 1 

.  0036 

6.2 

2330 

327 

129 

7.0 

12.3 

15.2 

21 . 3 

.  0035 

6.1 

0000 

329 

129 

6.8 

12.2 

15.2 

21.3 

.  0035 

6.0 

0030 

331 

123 

6.5 

12.2 

15.2 

21 . 5 

.0033 

6.0 

0100 

333 

127 

6 . 3 

12.1 

15.2 

21.5 

.0033 

5.9 

0130 

335 

127 

6.1 

12. 1 

15.2 

21 . 6 

.  0032 

5.9 

0200 

337 

126 

6 . 1 

12 . 0 

15.2 

21.6 

.0032 

5.8 

0230 

339 

125 

6.2 

12.0 

15.2 

21 . 7 

.  0031 

5.9 

0300 

341 

125 

6.2 

12.0 

15.2 

21.7 

.  003i 

5.9 

0330 

343 

125 

6.2 

ii  .  9 

15.2 

21 . 9 

.  0032 

5.9 

0400 

344 

125 

6.2 

11 .9 

15.2 

21.9 

.0  032 

6 . 0 

0430 

346 

125 

6.2 

11 . 9 

15.2 

22. 0 

.0033 

6. 0 

0500 

348 

125 

6.2 

11.9 

15.2 

22.0 

.0  033 

6 . 0 

TABLE  II 

MABL  Model  Predicted  Atmospheric  Mixed  Layer  Quantities- 

Case  1 . 


HOUR 

Zi 

Zlcl 

WIND 

Th 

0430 

256 

243 

0.0 

13 . 3 

OSOO 

252 

216 

8.0 

13.4 

0S30 

250 

174 

8.0 

13.3 

0600 

250 

130 

8.0 

13.  0 

0630 

252 

117 

8.0 

12.9 

07C0 

256 

123 

8.1 

12.9 

0730 

259 

128 

8.1 

12.9 

0800 

262 

133 

8.1 

12.9 

0830 

265 

137 

8.1 

12.9 

0900 

267 

140 

8.1 

12.9 

0930 

269 

143 

8.1 

13.0 

1000 

272 

146 

8.2 

13.0 

1030 

274 

149 

8.2 

13.0 

1100 

276 

1S1 

8.2 

13.  0 

1130 

277 

153 

6.6 

13.1 

1200 

279 

153 

5.0 

13.0 

1230 

280 

iS3 

5.2 

13.0 

1300 

281 

152 

5.3 

12.9 

1330 

282 

152 

5.5 

12.8 

1400 

283 

ISO 

5.6 

12.8 

1430 

284 

149 

S  .  8 

12.7 

1S00 

285 

147 

6.0 

12.7 

1530 

286 

145 

6.1 

12.6 

1600 

288 

141 

6.3 

12.6 

1630 

289 

137 

7.0 

12.5 

1700 

291 

133 

7.8 

12.5 

1730 

294 

130 

7.7 

12. S 

1800 

296 

126 

7.7 

12.4 

1830 

299 

123 

7.6 

12.4 

1900 

302 

122 

7.6 

12.4 

1930 

305 

122 

7.5 

12.3 

2000 

2030 

308 

3it 

123 

124 

7.9 

8.3 

12.3 

12.3 

2100 

314 

126 

8.1 

12.3 

2130 

3i7 

126 

7.9 

12.3 

2200 

320 

127 

7.7 

12.3 

2230 

323 

127 

7.5 

12.3 

2300 

325 

128 

7.3 

12.3 

2330 

328 

128 

7.0 

12.2 

0000 

0030 

330 

332 

127 

127 

6.8 

6 . 6 

12.2 
12 . 2 

0100 

335 

127 

6.3 

12.1 

0130 

337 

126 

6.1 

12. 1 

0200 

338 

125 

6 . 1 

12. 0 

0230 

340 

125 

6.2 

12. 0 

0300 

342 

125 

6.2 

11.9 

0330 

344 

125 

6.2 

11.9 

0400 

346 

125 

6.2 

11.9 

0430 

347 

125 

6.2 

11.9 

0500 

349 

125 

6.2 

11 . 8 

SST 

MLD 

TW 

Ze 

15.2 

0.0 

0.0000 

0.0 

IS. 2 

0  .  0 

0  .  0000 

7.2 

15.2 

0.0 

0 . 0000 

6.9 

15.2 

0.  0 

0.0000 

6.4 

15.2 

0.0 

0.0000 

6 . 0 

15.2 

0.0 

0.0000 

5.9 

15.2 

0.0 

0.0000 

6.0 

15.2 

0.  0 

0.0000 

6.1 

15.2 

0.0 

0.0000 

6.2 

15.2 

0.  0 

0. oooo 

6.2 

15.2 

0.0 

0 . 0000 

6.2 

15.2 

0  .  0 

0.0000 

6.2 

15.2 

0.0 

0 .0000 

6 . 3 

15.2 

0.  0 

0  .  0000 

6.3 

15.2 

0  .  0 

0 .0000 

6.3 

15.2 

0.  0 

0.0000 

5.8 

15.2 

0.0 

0.0000 

5 . 2 

15.2 

0.  0 

0.0000 

5.2 

15.2 

0.0 

0.0000 

5.3 

15.2 

0.0 

0.0000 

5.4 

15.2 

0.0 

0 . 0000 

5 . 5 

15.2 

0.0 

0.0000 

5.6 

15.2 

0.0 

0 . 0000 

5.6 

15.2 

0.  0 

0 .0000 

5.7 

15.2 

0.0 

0 . 0000 

5.8 

15.2 

0.0 

0.0000 

6.0 

15.2 

0.0 

0 .0000 

6.3 

15.2 

0.  0 

0.0000 

6.2 

15.2 

0 . 0 

0.0000 

6.2 

15.2 

0.0 

0.0000 

6 . 2 

15.2 

0.0 

0 .0000 

6.1 

15.2 

0.  0 

0.0000 

6.1 

15.2 

0 . 0 

0 . 0000 

6 . 3 

15.2 

0.  0 

0.0000 

6.4 

15.2 

0 . 0 

0 .0000 

6.4 

15 . 2 

0.0 

0.0000 

6 . 3 

15.2 

0  .  0 

0 .0000 

6.3 

15.2 

0.  0 

0.0000 

6.2 

15.2 

0 . 0 

0.0000 

6.2 

i5 .2 

0.  0 

0 . 0000 

6 . 1 

15.2 

0 . 0 

0.0000 

6 . 0 

15.2 

0.  0 

0.0000 

5.9 

15.2 

0.0 

0 .0000 

5.9 

15.2 

0.  0 

0.0000 

5 . 9 

15.2 

0 . 0 

0.0000 

5.9 

15.2 

0  .  0 

0.0000 

5.9 

15.2 

0 . 0 

0.0000 

6 . 0 

15.2 

0.  0 

0.0000 

6 . 0 

15.2 

0.0 

0 .0000 

6 . 0 

IS. 2 

0.  0 

0. 0000 

6. 0 
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Figure  14.  Selected  Ocean  Vertical  Temperature  Profiles- 
Case  1 . 
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Fiaure  15.  Selected  Ocean  Vertical  Profiles  of  *:he 

Horizontal  Wind  Driven  Currents-  Case  i.  Fast 
west  component  u  (solid  line)  North-south 
component  v  (dashed)  . 
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TABLE  III 


Coupled  Model  Predicted  Ionospheric  and  Oceanic  Mixed  Layer 

Quantities-  caso  i. 


HOUR 

0430 

0500 

0530 

0600 

0630 

0700 

0730 

0300 

0830 

0900 

0730 

1000 

1030 

1100 

1130 

1200 

1230 

1300 

1330 

1400 

1430 

1500 

1530 

1600 

1630 

1700 

1730 

1800 

1830 

1900 

1930 

2000 

2030 

2100 

2130 

2200 

2230 

2300 

2330 

0000 

0030 

0100 

0130 

0200 

0230 

0300 

0330 

0400 

0430 

0500 


Zi 

351 

350 

351 
353 
357 
362 
368 
373 
378 
383 
388 
393 
399 
403 
408 
413 
417 
421 
425 
430 
434 
438 
443 
448 
453 
458 
464 
470 
477 
483 
49  i 
498 
505 
512 
519 
526 
532 
539 
545 
552 
558 
564 
570 
576 
58i 
587 
592 
597 
60  1 
606 


Zlcl 

WIND 

Th 

SST 

316 

0.0 

12 . 0 

15. 0 

284 

6.4 

12.0 

15.0 

249 

6.4 

11.8 

14.9 

222 

6.4 

11.7 

14.9 

208 

6.4 

11.7 

14.9 

205 

6.5 

11.7 

14.9 

209 

6.7 

11 .7 

14.9 

215 

6.9 

11.8 

14.9 

222 

7.2 

11.9 

14.9 

229 

7.5 

12.0 

14.9 

235 

7.8 

12.1 

14.9 

242 

8.7 

12.2 

14.9 

248 

8.5 

12.3 

15.0 

254 

8.3 

12.4 

15.0 

26  0 

8.2 

12.5 

15.0 

264 

8.1 

12.6 

15.0 

268 

8.1 

12.7 

IS.  0 

272 

8.2 

12.7 

15.0 

275 

8.3 

12.8 

15.0 

277 

8.5 

12.3 

15.0 

278 

8.7 

i2 . 9 

15.0 

279 

8.9 

12.9 

15.0 

279 

9.0 

13.0 

15.0 

278 

9.2 

13.0 

15.0 

276 

9.4 

13 . 0 

i5 . 0 

274 

9.6 

13.0 

15 . 0 

271 

9.8 

13.0 

15.0 

268 

10  .  0 

13.0 

15.0 

267 

10.5 

12.9 

15.0 

267 

11.0 

12.9 

15.0 

268 

10.9 

13.0 

15.0 

271 

10.8 

13.0 

15 . 0 

273 

10.7 

13 . 0 

1  J  .  0 

276 

10.6 

13 . 0 

15.0 

278 

10.5 

13 . 0 

15.0 

231 

10.4 

13.  1 

15.0 

283 

10.4 

13 . 1 

15 . 0 

285 

10 .3 

13.  i 

15 . 0 

287 

10.2 

13.1 

14.9 

289 

10.1 

13. 1 

14.9 

291 

10  .  0 

13 . 1 

14 . 9 

293 

9  9 

13 . 1 

14.9 

295 

9.8 

13. 1 

14 . 9 

296 

9.5 

13 . 1 

14.9 

298 

9.3 

13 . 1 

14 . 9 

299 

9 . 0 

13.  1 

14.9 

300 

8.8 

13. 1 

14 . 9 

301 

8.6 

13.1 

14.9 

302 

9 . 0 

13 . 1 

14 . 9 

303 

9.4 

13.0 

14.9 

HLD 

0.0 

22.5 
22.7 

22.7 

22.8 
22.8 

?§:§ 

22.9 

22.9 

22.9 

22.9 
20.1 
20 . 1 
20.0 
20.0 
16.8 
16.8 
17.3 

17.3 
18.7 
18.7 

20.5 

20.5 

22.4 
22.4 
23.1 

23.1 
23.3 

23.3 

23.6 

23.6 

23.9 

23.9 

24.2 

24 . 2 

24.4 
24.4 

24.7 
24 . 7 

24.9 
24 . 9 
25. 1 
25 . 1 

25.3 

25.3 

25.4 

25.4 

25 . 5 

25.5 


TU 

0.0000 
0.0000 
.0057 
.  0057 
.0032 
.0032 
.0019 
.0019 
.0004 
.0004 
-.0014 
-. 0014 
-.0028 
-.0028 
-.0038 
-. 0038 
-.0044 

-  .  0044 
-.0  042 

-  .  0042 
-.0030 
-.0030 
-.0013 
-.0  013 

.0007 
.  0007 
.0024 
.0024 
.0  037 
.  0  037 
.  0044 
.  0044 
.0045 
.  0045 
.0044 
.  0044 
.0  043 
.  0043 
.0  043 
.  0043 
.0  042 
.  0042 
.0041 
.  0041 
.0041 
.  0  04i 
.0  039 
.0039 
.0  038 
.  0038 
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TABLE  IV 


MAEL  Model  Predicted  Atmospheric  Mixed  Layer  Quantities- 


HOUR 

Zi 

Zlcl 

WIND 

Th 

SST 

MLD 

TW 

Ze 

0430 

351 

316 

0.0 

12.0 

15.0 

0.0 

0.0000 

0  . 

0500 

350 

284 

6.4 

12.0 
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ure  21.  Selected  Ocean  Vertical  Temperature  Profiles- 
Case  2. 


Figure  22. 


Selected  Ocean  Vertical  Profiles  cf  the 
Horizontal  Kind  Driver.  Currents-  Case  2.  Eas¬ 
iest  component  u  (solid  line}  North-south 
component  v  (dashed)  . 
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Figure  23.  Selected  Chanqe  in  Dceai  Vertical  Salinity 
Profiles-  Case  2. 
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ire  26.  Time  Evolution  of  th  e  Sea-Surf  ace  Temperature-:.. 

Case  3.  Predicted  values  are  shown  by  rhe  sond 
line.  Observed  values  are  shown  by  the  dashed 
line. 
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Figure  27.  Time  Evolution  of  Ocean  Mixed  Layer  Depth  (MLD) 
Case  3. 
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ure  28.  Selected  Ocean  Vertical  Temperature  Profiles- 
Case  3. 


Figure  29. 


Selected  Ocean  Vertical  Profiles  of  the 
Horizontal  Wind  Driver.  Currents-  Case  3.  Eas- 
wesx  component  u  (solid  line)  North-south 
component  v  (dashed)  . 
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